Frozen metastable states in ordered systems of ultrafine magnetic particles 
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For studying the interplay of dipolar interaction and anisotropy energy in systems of ultrafine 
magnetic particles we consider simple cubic systems of magnetic dipoles with anisotropy axes point- 
ing into the ^-direction. Using Monte Carlo simulations we study the magnetic relaxation from 
several initial states. We show explicitely that, due to the combined influence of anisotropy energy 
and dipole interaction, magnetic chains are formed along the z-direction that organize themselves in 
frozen metastable domains of columnar antiferromagnetic order. We show that the domains depend 
explicitely on the history and relax only at extremely large time scales towards the ordered state. 
We consider this as an indication for the appearence of frozen metastable states also in real sytems, 
where the dipoles are located in a liquid-like fashion and the anisotropy axes point into random 
directions. 

PACS numbers: 75.75.-(-a, 75.50.Lk, 75.40.Cx 75.40.Mg, 



In the last decade, systems of ultrafine magnetic 
nanoparticles have received considerable interest, due 
both to their important technological applications 
(mainly in magnetic storage and recordings) and their 
rich and often unusual experimental behavior [l| . An im- 
portant scientific question concerns the magnetic struc- 
ture of the systems |3|- Several experiments on disor- 
dered magnetic materials present indications of a spin- 
glass phase I3, S y, a S tD| or of a random anisotropy 
system [alSl: while on the theoretical side, it is a mat- 
ter of controversy if, at large concentrations of nanopar- 
ticles, a spin-glass phase exists or not llCi llll Ha . Il3l| . 
While Monte Carlo simulations on ageing [lOj and mag- 
netic relaxation [ll| seem to favorize the spin-glass hy- 
pothesis, simulations of the zero-field cooling and field- 
coolin g su sceptibility showed no indication of a spin-glass 
phase [13 ■ In thi s paper , we use Monte-Carlo simula- 
tions (see e.g. [ij, [ij, [l^) to study the slow magnetic 
relaxation from a non-equilibrium situation into the fi- 
nal state in ordered systems of ultrafine particles |l7j . 
We find that already in this fully ordered arrangement 
of magnetic dipoles, the competition between anisotropy 
energy and dipole interaction is sufficient to produce frus- 
tration and metastable frozen states. We consider this as 
an indication that in the corresponding real systems of 
ultrafine magnetic particles, where additional frustration 
naturally arises as a consequence of the disorder, spin- 
glass phases may exist. 

We focus on perhaps the most basic model of mag- 
netic nanoparticles that (i) assumes a coherent magneti- 
zation rotation within the anisotropic particles, and (ii) 
takes into account the magnetic dipolar interaction be- 
tween them. Here, in order to get insight into the in- 
terplay of dipolar interaction and anisotropy energy, we 
further simplify this problem drastically by (iii) placing 
all magnetic particles onto the lattice points of a sim- 
ple cubic lattice and (iv) orienting all anisotropy axes 
into the z-direction. It is known that the simple cubic 



dipolar system possesses a columnar antiferromagnetic 
(CAF) groundstate [igl, where the magnetic moments 
are arranged in linear chains along an arbitrary direction 
and each chain is surrounded by chains aligned in the op- 
posite direction. In our case, the additional anisotropy 
energy favorizes the formation of chains along the z-axis, 
which are arranged antiferromagnetically in the xy-plane, 
similar to the case of a cubic Ising system with additional 
dipolar interaction [ijj. 

For describing the magnetic structure, we are thus led 
to introduce two order parameters O^ and Ot, where 0^ 
describes the magnetic order in each chain and Ot the an- 
tiferromagnetic order between neighboring chains. Here, 
we study the relaxation of Oi and Ot , (i) from the CAF 
state [Oi = 1, Ot = 1), and (ii) from general configura- 
tions where all dipoles are randomly oriented [Og — 0, 
Ot = 0). We find that in contrast to Og, the transversal 
order parameter Ot depends strongly on the initial state. 
Only for relaxation from the CAF state, the state with 
lowest energy is reached quite fast, while for the gen- 
eral case, the system gets frozen in some intermediate 
disordered state. This dependence on the initial condi- 
tions becomes more pronounced in the thermodynamic 
limit and results, at low temperatures, in complex frozen 
structures that consist of several domains. In each do- 
main, the chains are ordered in an antiferromagnetic way. 

For the numerical calculations, we use the same model 
as in Refs. |ll[ll2l |. where every particle i of constant vol- 
ume V is considered to be a single magnetic domain with 
all its atomic magnetic moments rotating coherently. 
This results in a constant absolute value |/Xi| = MgV 
of the total magnetic moment of each particle, where 
Ms is the saturation magnetization. The energy of each 
particle consists of two contributions: anisotropy en- 
ergy and dipolar interaction energy. We assume a tem- 

(i) 

perature independent uniaxial anisotropy energy E^^ = 
—KV{{fiini)/\^i\)'^, where K is the anisotropy constant 
and the unit vector n.^ denotes the easy directions. As 



usual, the energy of the magnetic dipolar interaction be- 
tween two particles i and j separated by r^j is given by 
{fliflj)/rf, - 3{fl,nj){^ijfij)/rlj. Adding up the 
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two energy contributions and summing over all particles 
we obtain the total energy 
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In the Monte Carlo simulations wc concentrate on sam- 
ples oi N — L'^ particles placed on a cubic lattice with 
periodic boundary conditions and L between 4 and 12, 
where all anisotropy axes n^ point into the positive z- 
direction. The unitless concentration c is defined as the 
ratio between the total volume NV occupied by the par- 
ticles and the volume Vg of the sample. Here, we focus on 
the large concentration c/cq « 0.64, where cq = 2K/M'^ 
is a dimensionless material-dependent constant (cq ~ 1.4 
for iron nitride nanoparticles |llll2(]j), but we also tested 
samples with c/cq ~ 0.2 and 0.4 that gave qualitatively 
the same results. The relaxation of the individual mag- 
netic moments per Monte Carlo step is simulated by 
the standard Metropolis algorithm, where the pii are de- 
scribed by their two spherical coordinates 9i and ipi . The 
interaction energies are calculated using the Ewald sum 
method with periodic boundary conditions in x, y and 
z-direction (1%. il4i] . To study the magnetic relaxation we 
start from a given initial state and determine as a func- 
tion of time t (number of Monte Carlo steps) for each 
particle i the angle 9i between the magnetic moment /i^ 
and the z-axis, from which wc obtain the relevant quan- 
tities. 

To quantify the relaxation process and the final state, 
we study the dependence of the order parameters Oi{t) 
and Ot{t) on t (number of Monte Carlo steps). While 
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jTij \) describes the order along the chains, 

{Sj Sj+s) describes the order perpendicular to 

the chains. Here, m): is the z-component of the mag- 
netic moment of the jth chain and Sj = 1 if all mag- 
netic moments of the jth chain point into the positive 
z-direction, Sj = —1 if they point into the negative z- 
direction, and Sj — otherwise. The index j -I- ^ denotes 
a nearest-neighbor chain of chain j and ( ) denotes the 
average over all j', 5 and over the number Nc of config- 
urations that ranges from N^ = 1000 for the smallest to 
Nc = 100 for the largest systems. 

In order to identify a temperature above which the 
ferromagnetic order inside the chains and the antiferro- 
magnetic order between the chains cannot be preserved, 
we first consider relaxation from the CAF state. Fig- 
ures n^a,b) show Oi{t) and Ot{t) for systems of iV = 6'^ 
(white symbols) and 10'^ (black symbols) particles as a 
function of t, for several values of the reduced tempera- 
ture f = kBT/{2KV), where 2KV is the height of the 
anisotropy barrier. One can see that Oi{t) and Ot{t) 




1 kJ/(2KV) 2 



1 k,T/(2KV) 2 



FIG. 1: The order parameters (a) Oe{t) and (b) Ot{t) as a 
function of t (number of Monte Carlo steps) are shown for 
the system sizes L = 6 and L = 10 when starting in the 
CAF state for several values of the reduced temperature T — 
kBT/{2KV) = 10 (circles), 1 (squares), 2/3 (diamonds), 1/2 
(triangles up), 1/2.5 (triangles left) and 1/10 (triangle right), 
where T is the temperature, ks the Boltzmann constant, K 
the anisotropy constant and V the particle volume. In (c,d) 
Oi and Ot are shown for fixed t = 10'' as a function of T. 



reach final plateau values 0| and O^ quite fast. Fig- 
ures nic,d) show these values as a function of T. The 
figures show that for small temperatures, 0| and OJ" are 
close to 1 and decay rapidly to much smaller values in 
a narrow temperature regime around T w 0.5. In both 
cases, the width of the transition regime shrinks with in- 
creasing system size. From the inflection point, we can 
identify the critical temperature Tc ~ 0.5 above which 
the longitudinal and transversal order breaks down. The 
figures show that both, the order along the chains and the 
order in the xy-plane disappears at the same transition 
point. We can also see in the figure that only as long as 
T is below Tc, the order parameters do not depend on the 
system size L. For T > Tc, Ot decreases with L, which 
reflects the fact that the chains have started to break. 

Next, we consider the relaxation from a random initial 
conflguration. We start with T — 1/10, a temperature 
well below Tc. Figure [5] shows that again both, Oe{t) and 
Ot{t) behave in a similar way. They first increase with 
t and then reach a plateau value at the same crossover 
time ti{L) that increases monotonically with L (see the 
arrows in Fig. [21). But while Oi reaches an L- independent 
plateau value (which agrees with the equilibrium value 
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FIG. 2: The order parameters (a) Oi{t) and (b) Ot{t) are 
plotted versus the time t (number of Monte Carlo steps) when 
starting in random configurations for the temperature T — 
1/10 and for L = 4 (circles), 6 (squares), 8 (diamonds), 10 
(triangles up) and L = 12 (triangles left). The arrows indicate 
the approximate crossover times ti for the system sizes (from 
left to right) L = 4 to L = 12. 




FIG. 3: The order parameters Ot{t) are plotted versus the 
time t (number of Monte Carlo steps) when starting in ran- 
dom configurations for several L from L = 4 to L = 12 (same 
symbols as in Fig. |5J and for the temperature T = 1/5. The 
second crossover time i2 is shown for L — 4. 



Ol for T = 1/10 of Figs. nja,c)), the plateau values of 
Ot continue to increase with the system size and seem to 
approach the corresponding value O^ from Figs. nib,d) 
only in the thermodynamic limit. 

Therefore, we arrive at the following picture for T < 
Tc'. Below ii, the system is disordered, while above ti, 
all magnetic moments are ordered in chains along the 
z-direction. But in contrast to the relaxation from the 
CAF state (Fig.QJ, the chains arc not (yet) ordered in an 
antifcrromagnetic way. Accordingly, the system has not 
yet reached its (equilibrium) state of minimum energy 
and must be considered as metastable. We therefore ex- 
pect that there exists a second crossover time ^2, which 
can be roughly identified as the inverse of the proba- 
bility P{L,T) that a chain flips into a more favorable 
direction. Accordingly, above t2 '~ P^^{L,T) chain flips 
will dominate the relaxation and this leads to a further 
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FIG. 4: Visualization of the antiferromagnetic order in the 
a;y-plane for systems at /3 = 10 when starting in a random 
configuration, (a-d) One L = 10 system at fixed t (number 
of Monte Carlo steps), t = 10^ 10^, 10* and t = 10^ (e-h) 
one L = 4 system after t = 20, 10^ 2 • 10^ and 10^. The 
complete chains are indicated by + or — signs, depending 
on the direction of the chain. Antiferromagnetic domains are 
shown in black and white. Domains of the same color fit to 
each other and are allowed to emerge, whereas domain walls 
exist between clusters of different colors. Sites, where chains 
have not yet been built are indicated by the grey color, (e) 
The average relative size Nmax/L,'^ of the largest cluster at t 
with ti < i < t2 is plotted versus the size L for T = 1/10 
(squares) and T — 1/20 (diamonds). 



increase of Ot{t) towards its equilibrium value O^ (see 
Fig.^d)). Since P{L, T) decays drastically with decreas- 
ing T and increasing L, one has to go to small system 
sizes and comparatively large temperatures (below Tc) 
in order to find t2- Indeed, Fig. 13 shows that at the in- 
termediate temperature T = 1/5, i2 can be observed only 
for the smallest system size L — A. We like to note that 



the metastable state shown here differs from the typical 
metastable spin-glass structure identified in [l9| for an 
Ising system, where the up- and down-columns are ar- 
ranged in a completely random fashion and Of should 
be equal to zero. Instead, we will show in the following 
that in our system, different domains grow from different 
seeds. 

We now consider the L-dependence of Ot{t). The in- 
crease of the plateau values with L in Figs. Efb) and El 
seems to suggest that in the thermodynamic limit, the 
system reaches the CAF groundstate already at ti{L). 
To see that this is not the case, we illustrate the relax- 
ation process in Fig.^for /3 = 10 and two different system 
sizes L = 4 and L — 10. To visualize the antiferromag- 
netic order, we follow the definition of Ot'- each of the 
L^ sites in the xy-plane can be either a "-I-" site or a "-" 
site, if all magnetic moments in the chain point into the 
positive or negative ^-direction, respectively, or a " 0" site 
if this is not yet the case (grey sites), in close analogy to 
the definition of Sj. The figure shows that in the initial 
time steps, for both sizes, only distant isolated chains 
have been formed. Due to the dipolar interaction, chains 
are more likely formed in the neighborhood of another 
chain. This way, small (columnar) antiferromagnetic do- 
mains develop that grow further with increasing time. 
When different domains contact each other, they either 
emerge and form larger domains (if they fit to each other) 
or establish quite stable domain walls between them. In 
Fig. ^ both types of domains are in black and white, re- 
spectively, such that domains of the same color fit to each 
other and are in principle allowed to emerge. The figure 
illustrates the change of the domains (i) with t (number 
of Monte Carlo steps) and (ii) with system size L. 

First, Fig.0]shows that below ii, the longitudinal order 
is not yet fully established and in rare cases (particularly 
for small system sizes) complete chains can still flip into 
the other direction (compare Figs.QJf) and (g)). Above 
ii, chain flips are too rare to be observed and the struc- 
ture is practically frozen-in, until at considerably larger 
time scales (above ^2, not shown here) chain flips may 
again occur. However, since ^2 increases drastically with 
L, the structure is practically frozen-in above ti for rea- 
sonably large system sizes. The reason for this glassy-like 
phenomenon is the competition between the longitudinal 
order along the chains and the transversal order in the 
plane that leads to frustration of the single magnetic mo- 
ments inside the chains. 

Second, also the monotonous increase of the plateau 
value of Ot with L in Figs. El^b) and 13 can be under- 
stood qualitatively from Fig. ^ According to the fig- 
ure, the system consists of large clusters and of small 
inclusions inside them. Above ti, the fraction of the do- 
main wall sites decreases with increasing L, and thus 
Ot is enhanced. (Compare Figs. 0fh) and^Id): in the 
first case, all sites except one belong to a domain wall, 
whereas in the case of L = 10 many interior sites ex- 



ist.) Nevertheless, the system does not reach an homoge- 
neous antiferromagnetic order. Instead, as it is shown in 
Fig. ^i), we have found numerically that at small tem- 
peratures, the fraction of sites occupied by the largest 
cluster, Smax/L^, has a size-independent value that is 
close to 0.6 for f = 1/10 and 1/20. This means that the 
relative size of the largest cluster does not grow with L 
on dispense of the smaller clusters, i.e., no unique domain 
is reached. 

From the preceeding it is evident that different initial 
conditions will lead to different domain structures in the 
plateau regime. When we start e.g. with a random chain- 
like structure consisting of certain (columnar) antiferro- 
magnetic domains, the first crossover time ti vanishes as 
for the relaxation from the CAF and the initial state will 
freeze in. Therefore, in the plateau regime, the struc- 
ture of the frozen configuration depends strongly on the 
history. 

In summary, we have shown that already an ordered 
system of ultrafine magnetic particles where each parti- 
cle is located at a lattice site and its anisotropy axes are 
oriented parallel to each other, shows complex dynamical 
behavior with the formation of frozen history-dependent 
states, where linear magnetic chains are formed that are 
quite stable and act as seeds for the formation of the 
frozen domain structure. Accordingly, the complex be- 
havior results from the interplay between dipole interac- 
tion and anisotropy energy. We expect that also in real 
sytems of ultrafine magnetic particles, where the dipoles 
are located in a liquid-like fashion and the anisotropy 
axes point into random directions, similar linear struc- 
tures are formed in the beginning of the process and act 
as seeds for the formation of larger frozen-in structures in 
the same way as the chains in the ordered systems. We 
thus take the results for the ordered magnetic structure 
as an indication that also in the corresponding disordered 
system, where frustration arises in a quite natural way, 
spinglass phases exist at low temperatures. 
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